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We investigate the relation of the semigroup probability density of an infinite activity Levy 
process to the corresponding Levy density. For subordinators, we provide three methods to 
compute the former from the latter. The first method is based on approximating compound 
Poisson distributions, the second method uses convolution integrals of the upper tail integral 
of the Levy measure and the third method uses the analytic continuation of the Levy density 
to a complex cone and contour integration. As a by-product, we investigate the smoothness of 
the semigroup density in time. Several concrete examples illustrate the three methods and our 
results. 
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1. Introduction 

For the infinitely divisible laws, there are a number of intriguing and useful relations and 
points of similarity between the probability measures or probability densities of the laws 
on the one hand and their associated Levy measures or Levy densities on the other. 

In particular, if U is the Levy measure of an infinitely divisible law on M. d with associ- 
ated Levy process {X t }t>o and if P(dx;t) denotes the law of X t , then (see Sato (1999), 
Corollary 8.9) 

limt- 1 f f(x)P(dx;t)= f f(x)U(dx) (1) 

for any function / in the space C# of bounded continuous functions on R d vanishing in 
a neighborhood of 0. It may also be noted that Burnaev (2006) gives two formulas for 
computing the Levy measure from the corresponding cumulant function. 

The present paper considers the opposite problem, that of calculating P(dx; t) from 

U. 
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In many concrete examples, there is a power series representation 



00 j.n 

P(dx;t) = y2-rU n (dx), (2) 



n=l 



where the U n (dx) are, in general, signed measures and U\ is equal to U , the Levy measure 
of P(dx;t). When d = 1 and both P(dx;t) and U(dx) are concentrated on the positive 
half- line R>o = (0, 00), it is convenient to give (2) the form 



00 ,. 
t 



P+( X ;t) = J2-U+(x), (3) 
n—l 

with P + (x;t) = P((x,+oo);t) and U+(x) = U„((x,+oo)) denoting the upper tail inte- 
grals. 

More particularly, when P(dx,t) (resp. U(dx)) is absolutely continuous with density 
p(x; t) resp. u(x) - the setting on which we focus in this paper - then the 'density version' 
of (3) is 

p(x;t)=^2—u„(x). (4) 
n=i n - 

In (4), necessarily, ui = u and the question is how the further coefficients u n may be 
calculated from u, possibly also using properties of the cumulant function of X\ (which, 
of course, is essentially determined by u). As part of the problem, we discuss conditions 
ensuring that an expansion of the form (4) exists. This issue is of some independent 
interest. 

Note that in the case of a finite Levy measure, the process X is a compound Poisson 
process and validity of formulae (4) and (2), with straightforward modifications for the 
atom at zero, is easily established; see Barndorff-Nielsen and Hubalck (2006), Section 4.1. 
We shall not consider this case here any further. 

Except for the discussion in Section 2 and some remarks in Section 4, we only consider 
the case where the process X is a subordinator with infinite Levy measure and without 
linear drift. In other words, X is an infinite activity pure jump subordinator. 

We shall discuss three methods for determining the coefficients u n {x). The first in- 
volves, as the final step, a limiting operation. We have u n (x) = lim e ^o u ne (x), where 

u ne (x) =J2(-l) n - k (f) c(eT- k u* E k (x). (5) 

k=l ^ ' 

Here, u e {x) is an approximation of the Levy density u(x) that corresponds to a com- 
pound Poisson process with intensity c(e) and *k indicates fc-fold convolution. The second 
method uses derivatives of convolutions, namely, 

JTl 

u n (x) = (-l) n — ((U+D(x). (6) 
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In words: we can obtain u„(x) as the nth derivative of the nth convolution power of the 
upper tail integral of the Levy density u(x). 

The third method uses the complex contour integral 



where n(9) is the analytic continuation of the cumulant function to a complex cone con- 
taining dt6 > and the contour is, roughly speaking, along the boundary of the cone. We 
will see that such an analytic continuation can be derived from an analytic continuation 
of the Levy density u{x) to a complex cone containing the positive real axis. 

We proceed to mention various works containing results that are related to those of the 
present paper. The comprehensive monograph by Sato (1999) contains many instances 
of the interesting relations between the probability distributions and the Levy mea- 
sures of infinitely divisible laws (cf. also Embrechts et al. (1979), Embrechts and Goldic 
(1981) and Sato and Steutel (1998)). Some examples are the relation between unimodal- 
ity properties of the two types of densities (Sato (1999), Section 52) and the behavior 
under exponential tilting (or Esscher transformation). See also Leandre (1987), Ishikawa 
(1994) and Picard (1997) who, partly in the wider setting of pure jump processes, study 
cases where the transition density exists and behaves as a power of t for t — > 0. Continuity 
of P + (x; t) at t = is characterized in Doney (2004). In Riischcndorf and Woerner (2002) 
(also see Woerner (2001)), the authors have established the validity of expansions for 
the probability density or distribution function of X t that are related to, but essentially 
different from, (3) and (4). 

The paper is organized as follows. Section 2 consists of a number of initial remarks on 
the problem at hand. Section 3 contains our main mathematical results. Illustrative ex- 
amples will be given in Section 4. Technical auxiliary material used in the proofs appears 
in the Appendix. 

The results in the present paper build partly on our previous, unpublished preliminary 
work (Barndorff-Niclsen (2000) and Hubalek (2002)). 

2. Initial remarks 

2.1. A first, motivating example: the positive cn-stable 
distribution 

We consider the positive a-stablc distribution with Levy density 




(7) 



.r 



-l-Q 



u(x) 



n- a y 



(8) 



where < a < 1. Note that we interpret T(s) 1 as an entire function with zeros at the 
non-positive integers. Alternatively, we could use the functional equations T(s + 1) = 
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sr(s) and T(s)r(l — s) = 7tcsc(7Ts) to rewrite expressions in a more familiar (and more 
lengthy) form. 

In general, there is no closed form expression for p(x;t) in terms of elementary func- 
tions, but it is well known (see, e.g., Feller (1971), XVII. 7 (6.8)) that 



u n (x) 



n>l 



(-1)" , 

r(— na) 



-l—na 



(9) 



Let us first illustrate the calculation of U2(x) by our first method, that is, as limit 
of U2 £ {x) for e — ► 0, where we use the approximation u e (x) = h e>00 \{x)u(x). We have 
c(e) = E~ a /Y{1 — a). For x > 2e, we obtain, by symmetry and partial integration, 



*2i \ 

u e (x) 



ar(-a) 



e~ a {x-e) 



-l-ct 



-1— 2a nx/2 

+ (l + a)/ y- a (x-y)- 2 - a dy 



(10) 



As U2 £ (a;) = u* 2 (x) — 2c(e)u e (x) , we obtain, in the limit, 



1*2(2) 



oT(-a) 2 



(1 + a) 



x/2 



v a (x-y) 



-2-a 



dy 



(11) 



The integral on the right-hand side can be expressed in terms of the incomplete beta 
function and, in this particular case, reduced to integrals for the complete beta function 
by elementary substitutions. This finally yields agreement with (9) for n — 2. 

In principle, though less explicit and more cumbersome, the method can be used 
for n > 2. Instead, let us illustrate our second method, the calculation of u n (x) according 
to formula (6). The tail integral is U + (x) — x~ a /T(l — a) and, by induction, or, quicker, 
by looking at the Laplace transforms, we see that (U + )* n (x) = x n ~ na /T(n — na). Dif- 
ferentiating this equation n times and applying the functional equation of the Gamma 
function to simplify the expression, we obtain (9). 

Finally, let us illustrate our third method, the calculation of u n (x) according to formula 
(7). The Laplace cumulant function is k(6) = —9 a and we get 



u n (x) 



(-1)" 
2m 



(12) 



To see that this, in fact, gives (9), we have to substitute 9 > 0/x and recognize the 
resulting integral as a variant of the Hankel contour integral for T(—na) . 



2.2. A simple, general result 



Let P(x; t) and u(x) be, respectively, the cumulative distribution function and the Levy 
density (assumed to exist) of an infinite activity Levy process on K>o- Let the u e (x) be 
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integrable Levy densities that we think of as approximations of u(x) and let us define 
c(e) and u ne {x) as in (5) above, setting 

/>oo 

U 0e (x) = l, U„ e (x)=- u ne (y)dy (n>l). (13) 



Theorem 1. Suppose 



Then 



lim / (1 Ax)\u £ (x) -u(x)\dx = 0. (14) 



t n 

P(x;t) = ]imJ2u ne (x)- i , (15) 

fr-^fl ^ < 17.1 



n>0 



pointwise for each x € M>o and t > 0. 



This follows easily from Levy's continuity theorem and the observation that the 
distribution of an infinite activity process is continuous. Details can be found in 
Barndorff-Nielscn (2000) and Barndorff-Nielsen and Hubalek (2006). Note that the ap- 
proximation u e {x)I[ e ^ ao ) satisfies the assumptions of Theorem 1. 

2.3. Formulae for u nE 

The first few functions u nE are 



Ul e {x) 


= Ue(x), 


U2e{x) 


= uf(x) 


U3e(x) 


= uf(x) 


U4 e (x) 


=< 4 -' 



[e)uf(x) + 3c(e)\ 

,*3 i R„fr\2„.*2f„\ 



(16) 



Further, it follows from the well-known inverse relations for binomial sums (sec, e.g., 
Comtet (1970), III. 6. a e) that the formula defining u n£ (x) in (5) implies conversely that 



uT(x) = J^{l)c(er- k u ke (x). (17) 

k=l ^ ' 

This can be used to compute u ne {x) inductively by 

n-i / v 

u n£ (x) = u?(x) -E(fc c(eT- k u ke (x). (18) 



fc=i 
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In particular, we have 



uu(x) 


= u e (x), 






u 2e {x) 


= uf(x)-2c(e)u u {x), 








= uf(x) - 3c(e)u 2e (x) - 


- 3c(e) 2 w £ 


(*). 


u 4e {x) 


= ii* 4 - 4c(e)u 3e - 6c(e] 


2 u 2 e(x) - 


Ac{efu e (x) 



(19) 



2.4. Cancellation of singularities 



Convergence of u ne (x) to a function u n (x) implies a subtle cancellation of singularities. 
Equation (5) contains an alternating sum of terms that diverge. Both c(e) and the con- 
volution powers of u £ {x) tend to +oo as e — > 0; in particular, see the first few instances 
of that formula listed in (16) above. 

To gain an understanding of how this cancellation occurs, note that for n > 0, we have 

Un+le (:r) = (n+l)a ; - 1 jY u n£ {x - y)u e (y)dy + {-l) n c{e) n u e {x)^, (20) 

with u £ (x) = xu e (x), as can be established by induction. Consider the case n = 2, and 
let U+(x) = J^°u e (y)dy and U + (x) = J^°u(y)dy. Using (20) and noting that c(e) = 
Jq u e (y) dy + (x) , we may rewrite u 2 e{x) as 



u 2e (x) = 2x jy u s (y){u e (x - y) -u e (x)}dy-u e {x)U+{x)\. (21) 

Hence, by letting e — > and invoking condition (14), we obtain the following. 

Proposition 2. Suppose the Levy density u(x) is differentiable and let u(x) =xu(x). 
Then 



u 2 (x) = 2x- 1 ^J^ u(y){u(x-y)-u(x)}dy-u(x)U + (x)^, (22) 

with the integral existing and being finite. 

Formula (21), first given in Barndorff-Nielsen (2000), has been generalized by Woerner 
(2001) to 



1 - n 1 

u n+ u{x) = - 



71+1 



u e (y){une(x -y)- u n£ (x)} dy - u n£ (x)U + (x) 



(23) 



where u ne (x) = xu ne (x). Typically, though, this cannot be used to pass to the limit e — ► 
for n > 1 . 
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2.5. Miscellaneous further points 

• The power scries representation (4) with the coefficients u n (x) given as limit 
of u n£ (x) was originally derived by a heuristic argument that may 
be found in Barndorff-Nielsen and Shcphard (2008), Section 3.3.3; sec also 
Barndorff-Niclscn and Hubalek (2006), Section 2.1. In Section 3 below, we give a 
rigorous derivation for the case of one-dimensional subordinators, under fairly strong 
assumptions. In subsequent work, we hope to establish proofs for Levy processes in 
R d . 

• Heavy tails do not matter much in the problems studied in the present context as we 
can deal with them using the Esschcr transform. Rather, it is the behavior of small 
jumps that can cause difficulties, as reflected in the assumptions of the theorems 
given below. For details, see Barndorff-Nielsen and Hubalek (2006), Section 2.4. 

• Suppose p(x; 0) = and 

limt~ 1 p(x;t)=u(x), (24) 

and let u £ (x) = e~ 1 p(x; e), in which case c(e) = e . By using the semigroup property 
p* k (x,e) =p(x,ke), formula (5) takes the form 

u ne {x) =e-f>l)"- fc ( n k )p(x;ke). (25) 
fe=i ^ ' 

Now, the right-hand side of (25) is, in fact, an nth order difference quotient of p(x\ t), 
so, provided p(x; t) is n times diffcrcntiable from the right at t = 0, we have 

gn 

\imu ne (x) = —p(x;0). (26) 

Thinking of u e (x) as an approximation of e p(x;e), as well as an approximation 
of u(x), this is then a further indication that, in considerable generality, univariate 
as well as multivariate, p{x; t) may be calculated via (5) in the manner discussed 
above. Of course, in practice, choosing u £ (x) = e p(x;e) is not an option since the 
point is to determine p(x;t) in terms of the Levy density u(x). 

• As a referee pointed out to us, our work has a strong connection to semigroup theory, 
in particular to the so-called exponential formulae. Let T be the transition operator 
defined by 

/>oo 

Tf(x)= f(x + y)p(y;t)dy. (27) 



For example, Theorem 1, which holds true for infinitely divisible distributions on 
R d , can be seen as a variant of Hillc and Phillips (1957) (£2), page 354. Yet, our 
main interest is in interchanging the delicate limit e — > and the infinite summation, 
and we have not been able to use results from semigroup theory for that purpose. 
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• Condition (14) is satisfied, in particular, if there exists an integrable function v on 
M>o such that (x A l)u e (x) < v(x) for all x € R>o and all e. Some candidates for u £ 
are 

u s (x) = lr £i00 )(x)u(x) or u s (x) = u(x)e~ x ^ £ . (28) 



3. Main results 

In this section, we analyze the following issues. 

• Does u ne (x) converge as e — > 0? If so, can wc find a more direct method to compute 
the limit u n (x) from it(x) and so avoid the difficult cancellations in u nE [x) as e — * 0? 

• If we have convergence, is p(x; in fact n-times diffcrentiablc (from the right) at 
t = and, if so, is w„(i) the nth derivative? 

• If the answer to the previous question is yes for all n > 1, do we have a convergent 
Taylor expansion of p(x;t) at t = 0? Is p(x;i) in fact an entire function in t £ C? 

In the proofs, we refer to several technical estimates that are provided in the form of 
lemmas in the Appendix. Recall that we are assuming that the process X is an infinite 
activity subordinator. 



3.1. Pointwise convergence of the coefficient functions 

In this subsection, we investigate the limiting behavior of u ns {x) as e — > for the partic- 
ular choice u £ (x) = c~ e / x u{x). This approximation is simple, it is always feasible, in the 
sense that it implies 



lim u £ (x) = u(x), lim / u e (y)dy~ / u(y)Ay Vx > 0, (29) 
e^O z^OJx J x 

and u e (x) will be smooth if u{x) is smooth, a property exploited below. We provide 
conditions on u(x) that imply the convergence of u ne (x) and obtain the expression (6) 
for the limit u n (x). 

Theorem 3. Suppose n G N and 

e- r V' +1 |u (fe) (2;)|dx<oo, fc = 0,l,... ! m ! (30) 



holds for some m>n + 2 and r > 0. 
(i) Let 

U+(x)= u(y)dy, x>0. (31) 

J X 

Then the nth convolution power of U + (x) is well defined for x > and (U + )* n (x) is 
n-times continuously differentiable. 
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(ii) If we set u s (x) = c~ £ / x u(x) for e > and x > and define 

u ne {x) ( n \ uf{x)e{sT- k , (32) 

fc=i ^ ' 

where c(e) = J °° u e (x) dx, then lim e ^o Une(x) = u n {x) for x > 0, where 

H™ 

u n {x) = {-lY — (U + T n {x). (33) 

Proof. In Lemma A. 9, we show that (U + )* n (x) exists, and in Lemma A. 10, that it is 
n-times differentiable. Smoothness of convolutions is less obvious than it at first seems 
(cf. Doctsch (1950), Section 2.14, page 104ff and Uludag (1998)). 

Let X ne {9) denote the Laplace transform of u ne (x). It is given by X ne (9) = K e (9) n — 
(— l) n c(e) n . Let X n (8) = n(8) n . Note that X n (8) is not the Laplace transform of u n (x); 
the Laplace transform of u n (x) does not exist. But, using the estimates from Lemma A. 5, 
namely 

1X^(0)1 <E mn /\6\ m - n , \X^\9)\<E mn /\Br- n (34) 

for some constants E mn , we see that both x m u ne {x) and x m u n (x) have integrable Laplace 
transforms, namely XhS (0) and A„ (#). According to Lemma A. 10, we can write the 
inversion integrals 

-i />r+ioo -i />r-j-ioo 

Un ^ x) ^2^J r _ ioc A ^W*d0, u n (x) = ^—^ /^ X^(°)e 6x de, (35) 



and obtain 



P + OO 

\u ne (x) - u n {x)\ < — / |A&>(r + iy) - X^(r + iy)| dy. (36) 



'ne 

oo 



We have < u s {x) < u(x), for x > 0, e > and u e (x) — > u(x) for e — > 0. Thus, looking at 

Ke {6)= / (c- 0x - l)u e (x)dx, k{6)= / (e - *" - l)u(x) dx (37) 
Jo Jo 

and 

/>oo /* oo 

4 k \9) = (-l) k e- 9x x k u e {x)dx, ^ k \9) = {-l) k e- ex x k u{x)dx (38) 
Jo Jo 

for k = l,...,n, we see, by dominated (or monotone) convergence, that lim 6 ^o K e($) — 
n(9), lim e ^ ^e*' (^) = K.( k \9). The functions A^(#) and A^, Tn '(#) are polynomials in 
K £ {9), K' £ (9),... 1 4 m) {9), respectively in k{9), k'{9), . . . , ^ m \9), thus \mx e ^ Q X^\6) = 
X^\8). Moreover, they are dominated by the integrable function E mn /\9\ m ~ n and by 
dominated convergence in (36), we have the desired result. □ 
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An interesting class of infinitely divisible distributions on R>o is the family of gener- 
alized gamma convolutions. This class is characterized by having absolutely continuous 
Levy measures with densities u(x) such that u{x) = xu{x) are completely monotone 
functions (Bondesson (1992), Theorem 3.1.1). 

Theorem 4. If u(x) is a Levy density such that u(x) = xu{x) is completely monotone, 
then the integrability assumptions (30) in Theorem 3 hold for all neN with arbitrary 
r>0. 

Proof. By the Bernstein- Widdcr representation for completely monotone functions, we 
know u(x) is holomorphic in 5ftx > 0, thus the Taylor scries at x > has radius of con- 
vergence x and 

oo n 

u(x/2) = u(x - x/2) = (- 1)" u {n) (x) ^ . (39) 

In our setting, J °° e~ 8x u(x/2) dec < oo for 8 > and as (— l)"^"' (x) > 0, we can integrate 
the series term by term. As u^ n \x) — xu^ n \x) + nu'" _1 '(i) for n > 1 and we know that 
J e~ rx xu(x)dx < oo, we inductively obtain the result. □ 

Remark 5. An example where Theorem 3 applies but the corresponding distribution 
is not a generalized gamma convolution is given by u(x) = x ~ 3 / 2 c sln ( x \ An example 
where the integrability conditions (30) do not hold for n > 1 and any r > is given 
by u(x) = x~ 3 / 2 sin(x~ 3 ) 2 . We do not know whether the conclusion of the theorem is 
nevertheless true in this case. 

3.2. Differentiability in time 

For the proof of differentiability properties of the probability densities p(x; t) with respect 
to t > 0, we need slightly different integrability properties of the cumulant function k(9) 
and its derivatives. Sufficient conditions to guarantee those from assumptions on the Levy 
density u(x) are conveniently formulated in terms of the integral modulus of continuity. 
We recall, for example, from Kawata (1972), Theorem 2.7.4, that the integral modulus 
of continuity uj^(S; f) for an integrable function f(x) and a real number 5 > is defined 

by 

/+oo 
\f(x + h)-f(x)\dx. (40) 
-oo 

We note that it is sufficient to consider < h < 5 in (40). We use the integral modulus of 
continuity for functions f(x) that are a priori defined for x > with the understanding 
that f(x) = 0, if x < 0. 



774 O.E. Barndorff- Niels en and F. Hubalek 

Theorem 6. Suppose meN, n £ N, a£ (0, 1), r £ [0, oo) and u(x) is the Levy density 
of an infinite activity subordinator. Suppose 

1 + na . . 

m>- , 41 

1 — a 

u(x) is m-times differentiate in x > 0, the functions 

v £ {x) = (-l) e c~ rx x e+1 u^(x) (£ = 0,...,m) (42) 

are integrable and their integral modulus of continuity satisfies 

cjM(5;vz) = 0(5 1 - a ) (<y-»0). (43) 

Let p(x;t) denote the probability densities corresponding to u(x). Then p(x]t) is, for all 
x > 0, n-times differ entiable in t> 0; furthermore, 

«*(*) = (- 1 ) fc ^(C /+ )* fc W (fc = 1, • • ,»), (44) 



zs we// defined and 



;i ,(x,0)=u k (x) (fc = l,...,n). (45) 



Proof. Let 



A " m)( ' ;<) = 9^ GK(r P " (X;<) = (46) 

We will show the following statement inductively for n' = 1, . . . ,n. We have, for all x > 0, 
that p(a;;f) is n'-times diffcrentiable in t > and that 

C — IV" r c + ioc , . 

p "' (x; ^ )= ki^y c _ ioo ^v-d* ( 4? ) 

First, with n' = 0, we observe that (— l) m 2: m p(x; i) has Laplace transform X^ n \9;t). In 
Lemma A. 8 below, we show that the assumptions on the integral modulus of continuity 
imply \ { m) (6;t) = 0{\e\~ m ^-^) as 3(6>) -> ±oo. Since (41) implies m > 1/(1 - a), we 

have that Xq (9;t) is integrable and we can apply the inversion formula. This is all that 
was to be shown for n' — 0. Suppose, now, we have shown the claim for some n' — 1 and 
want to show it for n' . We can write 

h~ 1 (p n '-i{x;t + h) —p n t-\{x;t)) 

(48) 

( _-\\m rc+\oo 

= M^J h-^x^t + h^-x^m^de. 
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In Lemma A. 8 below, we show that the assumptions on the integral modulus of continuity 
imply h- 1 (X i jl 1 (e]t + h)-X^}_ 1 {e;t))=0(\e\ a ^ m+n '^- m ) as 3(0) -> ±oo. Now, (41) 
implies that the integrand in (48) is dominated by an integrable function and we can 
apply dominated convergence as h — ► 0. This shows that p n i-i(x;t) is differentiable with 
respect to t and its derivative is given by (47). This finishes the induction. To complete 
the proof, we observe that A n (0) = «(0)" = ^r[c K W] t=0 = a£ O) (0;O) and, in view of (46), 
(47) and (35), we indeed obtain p n (x;0) = u„(x). □ 

Remark 7. If we consider n = 1, then Theorem 6 provides sufficient conditions for the 
(pointwise) validity of 

u{x) = limt~V(^;*) ( 49 ) 

for x > 0. Our assumptions are quite different from those given by Woerner (2001), 
Corollary 2.3. 



3.3. Power series representation in time 

The purpose of this section is to show that, subject to some regularity conditions, the 
probability densities p(x; t) are analytic functions in t that can be represented by a power 
series of the form (4). To be able to do so, we assume that the Levy density u(x) is an 
analytic function satisfying some growth condition. 

Theorem 8. Suppose 

a>0, 0<a<l, /3>-l, 7>0, < ^ < ^ (50) 
and the Levy density u(z) is an analytic function in a domain containing 

VK = {zeC:z^0,|arg(z)| < tjj}. (51) 

Assume, moreover, that 

u(z)^az- 1 - a + 0(\zf) asz^OinW (52) 

and 

u{z) = e>(e TSRz ) nsz-toomf. (53) 
Then the cumulant function 

/>oo 

k(9)= / (e~ ex -l)u(x)dx (54) 
Jo 

admits an analytic continuation from {9 £ C : > 7} to {8 € C : 9 7^ 7, | arg(0 — 7)| < 
w + ip} that goes uniformly to as 8 — > 00 in {0 G C: ^ 7, | arg(0 — c)| < ^ + where 
c > 7 is arbitrary, but fixed. 
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Furthermore, p(x; t) is, for all x > 0, an entire function intdC and we have the power 
series expansion 

P(x;t) = V\ n (:r)— , (55) 



where 



Unix) = ~~~~7 / K{6) n e ax d9, (56) 



c 



2m 

with C the contour | arg(# — c)| = ip, 9 — being passed on the left. 



Proof. Let v(z) ~u{z) — az 1 a and X(9) = /™e 9x v{x)dx. Using 

(e~ ex -l)x- 1 ~ a Ax = T{-a)-9 a , (57) 



o 



we have 



k(9) = aT(-a) ■ 9 a + X(0)-X(O), (58) 



valid for W > 0. 

If we let V-(x) = v(e~ 1 ^'x)e~ i ^ , then we have the growth estimates V-(x) = C(|a;|' 3 ) as 
x — > and = O(exp((7cos^)x)) as oo. Thus, the Laplace transform A_(0) = 

J Q e~ V-(x) dx is absolutely convergent for $19 > jcosip and A_(#) — » uniformly as 
9 — > oo in 9f?# > ccos?/;; see Doetsch (1950), Satz 4, page 142 and Satz 7, page 171. 

Next, we show that X(9) = X-(9e~ 1 '^) for real 9 > 7. Suppose n, > 1 and let us integrate 
c~ 0z v{z) over the closed contour consisting of a straight line from n~ x to n, a circular arc 
from n to ne -1 ^, a straight line from ne~^ to n^e -1 ^ and a circular arc from n~ l e~ l1/ ' 
to n^ 1 (sec Figure 1). By Cauchy's theorem, this integral is zero. The estimates (52) and 
(53) show that the contributions from the circular arcs vanish as n — > 00 and we obtain 



X(9)= I e- 9x v(x)dx= / e- 9z v(z)dz (59) 

Jo Jo 



c 



^(e-^e^dx^A-^c-^). (60) 



A similar argument shows that the function v + (x) = v(e"^a;)e 1 * has Laplace transform 
A + (#) = L e~ 9x v + (x) dx, which is absolutely convergent for 5ft# > jcosip and satisfies 
X + (9) — ► uniformly as — > 00 in 3J0 > ccosip, and A(0) = X + (9c^') for real > 7. 

Looking at (58) reveals that the Laplace transform of p{x; t), namely e K ^ 1 , is integrablc 
on the vertical line = c. Thus, we can use the inversion integral 

-1 />c+ioo 

p(x;t) = — e K W t+fte dfl. (61) 

2™ Jc-ioo 
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Figure 1. Integration contour for the analytic continuation of n(0). 




Figure 2. The contour used to derive (62) and the final contour C. 



Let n > 1 and consider the integrand Q K ( ) t + 0x on the closed contour consisting of the 
vertical line connecting c — i • n and c + i • n, the circular arc with center c and radius 
n going from c + i • n to c + c^ n / 2+ ^ ■ n, the straight lines connecting c + e'^ 71 / 2 ^) • n, 
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c and c + e~^ n ' 2+ ^ ■ n and, finally, the circular arc from c + e~^ n ' 2+ ^ • n to c — i • n, 
see (Figure 2). By Cauchy's theorem, the integral is zero. Again, looking at (58) and 
the properties of the analytical continuation \(9) reveals that the integrand vanishes 
uniformly on the circular arcs as n — > oo and, by Jordan's lemma, we conclude that 

P (x;t) = — [ e^ e)t+ex d9. (62) 
2m J c 

On C, the linear term 9x dominates n(9)t as 9 — > oo. Consequently, (62) makes sense for 
any t £ C, in contrast to (61), where t > is required for convergence. We observe that 
taking t = yields 



1 

27ti 



e 6x d0 = 0; (63) 

c 



thus, p(x; 0) = 0, according to our convention above. We can differentiate (62) under the 
integral. Let us consider h £ C with \h\ < 1 and t £ C arbitrary. Using (62), we can write 
the complex difference quotient 

p(x;t + h)-p(x;t) 1 f c K ^ e)h - 1 



2m 



- .c K W t+8x dO. (64) 
7c " 



Again invoking the asymptotic behavior of k(9) as 9 — > oo on C, we can, by dominated 
convergence, prove the existence of the complex derivative dp(x;t)/dt for all f € C and 
the formula 

^ t p(x;t) = ±J c K(e)-e«W t+ex de. (65) 
It follows that p(x; t) is an entire function in t and 

<} = 2^i /- K{er ' "^^^^ d ^ (66) 



□ 



4. Examples 

4.1. The positive a-stable law 

The results for the positive a-stable law with < a < 1 have already been discussed in 
an informal way in Section 2.1. In view of the simple form of the Levy density (8), it 
is straightforward to check that the assumptions for Theorems 3, 6 and 8 are satisfied. 
Consequently, the results in Section 2.1 can be derived rigorously by our three methods. 
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4.2. The gamma distribution 

Suppose X is the gamma process, for which X\ has the law r(i/, a) with parameters 
v = 1 and a = 1 . The probability and Levy densities are 

p(x;t) = f ^x t - 1 e- x , u(x) = x- 1 e- x . (67) 

To illustrate our results, we choose the approximation u e (x) = x e u{x). We note that 
u e (x) = T{e)p{x,e) and thus 

n 

i-k I n 



x) = r(l + eT ■ s- n £(-!)"-* ? P (x; ke). (68) 



fc=i 

So, we are basically in the situation discussed in the third remark of Section 2.5: the con- 
vergence of u ne (x) to u n (x) is equivalent to the convergence of the nth order difference 
quotient of (67) at t = to the nth derivative from the right. The gamma probability den- 
sity is, for any x > 0, an entire function in t and the coefficients in the series expansion (4) 
are given by 

n-l / 

n I ;_i„ i_n— fc— 1 x 



u n (x) =x l e , klc k In 



k=0 



The numbers Ca, arise in the expansion T(l + z) 1 — J2 n >o c nZ n - They can be expressed 
explicitly as 

F n _i( 7l -C(2),2C(3), . . . , (-l)"- 2 (n - 2)!C(n - 1)) 



(n-l) 



with Y n the complete exponential Bell polynomials as given in Comtet (1970), III.3.C, 7 
the Euler-Mascheroni constant and £ the Riemann zeta function. 

Let us now imagine that we did not know (67). The function xu(x) is obviously com- 
pletely monotone and Theorem 4, and so Theorem 3, applies. Let us illustrate the calcu- 
lation of U2(x) by formula (33). The tail integral is U + (x) = E\{x), where E\(x) denotes 
the exponential integral; see Abramowitz and Stegun (1992), 5.1.1, page 228. A direct 
calculation of (U + )* 2 (x) is not very explicit. Let us write 

U + (x)=V(x)-L{x), L(x) = ]nx, V(x) = In a; + E x (x). (69) 

This decomposition is useful because L(x) is simple, while V(x) and its derivatives are 
intcgrable at zero. We have V(0) — —7 and V'(0) = 1. Equation (69) implies 

(U + )* 2 (x) = L* 2 {x) - 2(L * V)(x) + V* 2 {x). (70) 

Next, we observe L* 2 (x) = (In 2 x - 21na; + 2 - n 2 /6)x and thus [L* 2 (x)}" = (2]nx)/x. 
To compute the second derivatives of (L * V)(x) and V* 2 (x), we can interchange differ- 
entiation and convolution by the usual formulas; see Doetsch (1950), 2.14.5, page 115ff. 
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Namely, we use 

[V* 2 (x)}" = (V')* 2 (x) + 2V(0)V'(x) (71) 

and 

[(V * L){x)]" = (V" * L){x) + V'{0)L(x) + V(0)L'(x). (72) 

The convolution integrals on the right-hand side of (71) and (72) can be computed in 
terms of the exponential integral. Combining the three contributions, those terms cancel 
and we obtain 

u 2 {x) =2ar 1 e _x (lnx + 7), (73) 

in agreement with (68) above. 

Finally, what can we say about Theorem 8 in this case? In its present form, it does not 
apply since (52) is not satisfied, though formula (56) is correct. The cumulant function 
is k{6) = - ln(l + 9) and 

«n(aO = / ln(l + 6) n c 9x d9. (74) 

Agreement of this formula with (68) can be established by referring to the Hankel contour 
integrals for the derivatives of T(z) -1 at z= 1. 

4.3. The inverse Gaussian distribution 

The inverse Gaussian distribution IG{5 1 7) with 5=1 and 7=1 has a Levy density of 
the form 

u{x) = -^x- 3/2 e- x ' 2 (75) 

and the probability density is 

p ( x . t ) = _L te t a .-3/2 e -(t» a -»+*)/2 < (?6) 



Using the generating function for the Hermite polynomials H n (x) , namely, 

e™- t2 =J2H n (x)- r (77) 

n>Q 



we find 

TJ 

(78) 



i n (x) = ^-^x-^c-^H^ (| 
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Let us choose the approximation u s {x) = e~ £2 ^ 2x ^u(x). We recognize that this is a mul- 
tiple of p(x, e) and, again, showing the convergence of u ne (x) to u n (x) essentially reduces 
to a study of the nth order difference quotient of p(x;t) at t = 0. 

Let us look at the second approach, based on the tail integral. Again, xu(x) is com- 
pletely monotone. We have 

where erfc(:r) is the complementary error function; see Abramowitz and Stegun (1992), 
7.1.2, page 297. Let us illustrate the computation of u%(x). By looking at Laplace trans- 
forms, wc establish (U + )* 3 (x) = 2^c~ x / 2 {2 + x) - (2x 2 + 6x) crfc( 1 /f ) and differenti- 
ating — (U + )* 3 (x) three times, we obtain 

u 3 (x) = -^=x- 5 / 2 (x~l)c- x/2 , (80) 

\/27T 

in agreement with (78) above. Finally, Theorem 8 applies, the cumulant function is 
k{6) = 1 - a/1 + 26 and wc get 

Unix) = / (1 - vT+20) V 1 d9. (81) 

2m J c 

Agreement of this formula with (78) can be established as follows. First, we substitute 
^ i — ^ (f? — l)/2 and, expanding the integrand by the binomial theorem, we obtain a sum 
of Hankel integrals of the form (12) with a = 1/2, producing a sum of powers of x. Using 
the well-known explicit form of the coefficients of the Hcrmite polynomials shows (78). 



4.4. Examples on R: The Meixner and the normal inverse 
Gaussian distribution 

The following example is not covered by the standing assumptions in this paper, as 
the Meixner distribution (see Schoutens (2003), Section 5.3.10, page 62), is an infinitely 
divisible distribution on R, not on R>o- 

Let us consider the Meixner distribution with parameters /i = 0, 5 = 1. a = 1 and 
13 — 0. It has the density 

p(x;t) = ^ T{t + [x)Tit ~ [x) 



ix T{2t) 
This expression can be expanded in a series, 

p(x;t) = Y]u n (x)- { , \t\<\x\ 

z — ' IV. 

n>l 
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with 



where 



and 



Un{x) = xsinh( TO ) y - l(ai(x) ' ' ■ ' ' a - l(2;)) ' 



ai{x) = ip(ix) +i/j(-ix) + 21n2 + 27 



a n (x) = ^ (n) (ix) + - (-l) n 2"(n - l)!C(n) (n > 2). 

Here, Y again denotes complete exponential Bell polynomials, 7 is the Euler-Mascheroni 
constant, tp is the digamma function and £ the Riemann zeta function. Note, however, 
that here, p(x; t) is not an entire function in t due to the poles of the Gamma function. 
Thus, we must expect qualitative differences to the cases studied in the present paper. 

The normal inverse Gaussian Levy process with parameters /i = 0, 8 = 1, a = 1 and 
(3 = has the probability density 



p(x;t) = — U X— , - ieB. (82) 

This expression admits a power series expansion in t — for |t| < |x| and the coefficients 
u n {x) can be expressed in terms of polynomials and the Besscl K- function. We do not 
give details here, but note that this result is similar to the Meixner case. 

As both the normal inverse Gaussian and the Meixner Levy processes can be expressed 
as subordinated Brownian motions, it might be interesting to investigate the power series 
expansion of a subordinator and the corresponding subordinated process in general. This 
issue is left to future research. 



4.5. A bivariate example: The inverse Gaussian— normal inverse 
Gaussian law 

Again, this example is not covered by the assumptions of the rest of this paper as it deals 
with a bivariate distribution on K>o X R. We consider the probability densities 



P(x,y,t) = 2 exp 



r 



y 



(83) 



on R>o x R. They correspond to an inverse Gaussian-normal inverse Gaussian, or IG- 
NIG, Levy process. For properties of this type of law and its origin in a first-passage 
time problem for a bivariate Brownian motion, see Barndorff-Nielsen and Blaisild (1983), 
Example 4.1 and Barndorff-Nielsen and Shcphard (2001), Example 4.3. 
The associated Levy measure is 



u{x,y) 



1 

2n 



x exp 



2 I x 



(84) 
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and the Laplace cumulant function is k(9,t)) = 1 — y/l + 20 — rf. We can set u e (x,y) = 
u(x,y)c~ e /( 2a: ' and proceed, just as in Section 4.3 for the IG distribution, to obtain, 
with (a bivariate extension of) our first method, 

u n (x,y) = limJ2 (f) (-l) n - k e n - k u?{x,y). (85) 

£ ^ k=l ^ ' 

The resulting series expansion 
p(x,y;t) = ^ n (x,y), u n (x,y) = rtf-*"- 1 )/ 3 *"'"" 1 )' 2 ^. (Jfj (86) 

n>l U ' V V / 

agrees with the series obtained directly from the explicit expression for p(x, y\ t) and the 
generating function for the Hermite polynomials; see (77). 

Appendix: Auxiliary results 

This section provides technical estimates used in the proofs of Theorem 3 and Theorem 6. 
Several statements look well known and standard. Yet, for rigorous proofs, a careful 
checking of integrability and differentiability conditions is necessary and, in a few places, 
delicate, in the present setting. 

A.l. Auxiliary estimates for the cumulant function 

The structure of this subsection is as follows. In Lemma A. 2, we show that the integrabil- 
ity assumptions (30) of Theorem 3 imply a certain asymptotic behavior of the derivatives 
of u(x) as x — > and as x — ► oo. This is used in Lemma A. 3 to derive estimates of the 
derivatives of the cumulant function n(9) as 9 — > oo, by partial integration. We then 
introduce 

MO) = B nk ( K '(6), . . . , K ("- fc+1 >(0)) (87) 

and 

P* nk (6)=B nk (\K'(6)\,...,\K( n - k+1 H9)\), (88) 

where B nk denotes the partial Bell polynomials, as defined and discussed in Comtet 
(1970), Section 3, page 144ff. The estimates for the derivatives of k(9) are plugged into 
the Bell polynomials (3 n k{9) in Lemma A. 4. Using the latter, we obtain estimates for 
the derivatives of X n (9) as 9 — > oo in Lemma A. 5. Next, Lemma A. 6 shows that u E (x) 
satisfies the assumptions (30) uniformly for < e < 1 and thus, applying Lemmas A.2- 
A.5 to u £ (x), gives uniform estimates for the derivatives of A ne (#) as 9 — > oo. Finally, 
Lemma A. 7 provides a refined estimate for (3 n k{9) from the slightly stronger assumptions 
of Theorem 6. 
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We are considering distributions on K>o, thus we have |e K W| < 1 for 3?(0) > 0. More- 
over, k{9) is analytic for 9?(0) > and 

/>oo 

K {n) (9) = {-l) n C - ex x n u(x)dx, »(0)>O,n>l. (89) 
./o 



Definition A.l. Suppose n £ N and c > 0. FKe i/ien say i/iat assumption A n {c) holds 
for u{x) if u{x) is n-times continuously differentiate and 

e- cx x k+1 \u {k) (x)\dx<oo, k = 0,l,..,,n. (90) 

A consequence of assumption A n (c) is, that the Laplace transform J °° e~ 6x x m u^ n > (x) dx 
exists for any m > n + 1 and 9?(0) > c. 

Lemma A. 2. Suppose n>2 and c > 0. J/ assumption A n (c) holds, then 



lim e - te a;"u ( ™- 2) (s) = 0, lim e^cc V™~ 2) (x)=0 (91) 



/or 9?(0) > c. 



Proof. Let < a < b. Partial integration gives 

b r b 

e- ex {x n u {n - 2 \x))' dx^c- ex x n u^- 2 \x)\ b a +6 / c- ex x n u^- 2 \x)dx. (92) 



We have {x n u^ n ~ 2 \x))' = nx n ~ 1 u^ n ~ 2 '(x) + i™!i' n ~ 1 ' (x) and see from the integrability 
assumptions A n (c) and 3J(0) > c that both integrals in (92) converge to a finite value 
as a — > and 6 — > oo (separately). Thus, the limits lim^^o e~ ex x n u^ n ~ 2 ^ (x) =a and 
lim^^oo e~ 8x x n u^ n ~ 2 ' i (x) — to exist with finite a and w. But, a ^ or w ^ would imply 
that e _6 ' :r x™~ 1 u(™ _2 )(a;) is asymptoically equivalent to a/x as x — * 0, respectively to w/x 
asu oo. Both properties would contradict the integrability of e~ 9x x n ~ 1 u^ n ~ 2 \x) that 
follows, again, from assumption A n {c). Thus, we must have a = and ui = 0. □ 

The following lemma is essentially a reformulation of the well-known fact that an n- 
timcs diffcrcntiablc function f(x) with f^(x) integrable for < k < n has a Fourier 
transform f(y) which satisfies f(y) = 0(\y\~ n ) as \y\ — * oo. As we will need uniform 
growth estimates later, we provide a more detailed statement with explicit bounds. 

Lemma A. 3. Suppose n > and c > 0. If assumption A n (c) holds and we let 

/•oc 

L k (c) = / e- cx x k+1 \u^(x)\dx, fc = 0,l,...,n, (93) 



o 
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then 

|« (n) (0)|<^# m>c, (94) 

where 

M (c)=L (c), M n (c) = J2( n k ) (n)n-i-fci*(c) (n>l). (95) 

fc=0 ^ ' 

Proof. For n = 1, we have k'{9) = — J °° e~ ex xu(x) dx and the assertion of the lemma is 
obvious, namely, < Lq. For n = 0, we can write k(9) — J® k'(Q d£ and the assertion 

follows, namely |«(0)| < L o \0\. For n > 2, we recall k^(6) = (-1)" J °° e~ ex x n u(x) dx. 
Let < a < b. Repeated partial integration gives 

e- f)x x n u{x)dx^-^-^e- ex {x n u{x)) {k - 1) +^Zt/ Q~ 0X {x n u{x)) {n - 1) dx (96) 
fe=i a ^ a 



and, by the Leibniz rule, we obtain 

fc-i 



{x n u{x))^=Y,( k ^ 1 ) (nh-i-^-'-^WHx). 



(97) 



From assumption A n {c) and Lemma A. 2, we conclude, letting a — > and 6 — > oo, that 

K (™)(0) = i_L. y e-^Ca^ufs))^- 1 ^. (98) 
Using (97), this time with k = n, we get 

(x n u(x))( n -V = Y^ ( n ~ t 1 ) (n)n-i-tx e+1 uW(x). (99) 

This shows that the integral in (98) is bounded by M n {c). □ 
Lemma A. 4. Suppose n > and c > 0. J/ assumption A n {c) holds, then 

Pnk(0)<r^, k = l,...,n, (100) 

where M n k = B n f.(Mi,...,M n -k+i)> the constants M\,...,M n are as in Lemma A. 3 
above, and B n k denote the partial Bell polynomials. 
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Proof. The Bell polynomials B nk have non-negative coefficients and are therefore in- 
creasing functions of each argument. Using the bounds from Lemma A. 3, we have 

(3* nk {6) = B nk (\ K >(9)\,. . . , |) < B nk U,^,. . . , • ( 101 ) 

Using the homogeneity property of the Bell polynomials (sec Comtet (1970), Theo- 
rem HI. 3. A), in particular, the last part of the proof, we obtain the desired result. □ 

Lemma A. 5. Suppose m > 1, n > 1 and assumption A n (c) holds for u(x) with some 
c> 0, and let A„(0) = re(0) re . Then 

mAn 

\ X( ™H8)\ < E mn = J2 (n)jLr j M mj . (102) 

Proof. From the explicit form of Faa di Bruno's formula (see, e.g., Gradshteyn and Ryzhik 
(2000), (0.43) and Comtet (1970) (Theorems III.3.A and III.4.A)), we get A,^™' (0) = 
Y^jLx l ( n )j K (^) n ~'' x Pmj(O) and, in conjunction with the estimates from Lemma A. 4 we 
obtain the result. □ 

Lemma A. 6. Suppose n > 0, assumption A n (c) holds for u{x) with some c > and we 
let 

L k (c)= / e- cx x k+1 \u^{x)\dx, k = 0, 1, . . . ,n. (103) 



If we set u e {x) = c 6 / x u(x) for x > and e > 0, then assumption A n (c) holds for u e (x) 
and we have, for any e > 0, the uniform bound 

e- cx x n+1 \4 n] (x)\ dx < L n (c), (104) 



o 

where 

Z «( C )=EE(fc) {\Z])j/^ L n-u{c). (105) 
k=0t=l \ S \ / ' 

Proof. This follows from 

[e-^]( fc )=^(-l)^Y^: 1 )^f£)V^ (106) 



i=i 



e-ij tv 

and the inequality < x l c~ x < £ e e~ e for x > 0. □ 
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Lemma A. 7. Suppose m € N, aG (0,1), c £ (0, oo) and u(x) is the Levy density of an 
infinite activity subordinator which is m-times differentiable and such that the functions 

v e (x)= C - cx x e+1 u w (x) (£ = 0,...,m) (107) 

are integrable and their integral modulus of continuity satisfies 

J 1 \S;v l ) = 0(5 1 - a ) (6->0). (108) 

We then have, for n = 0, . . . , m, 

K ( n) (0)=O(\9\ a - n ), (3(0) — ±oo) (109) 

and, for I — 1, . . . , m, 

f3 me (6)=(D(\6\ ea - m ) (9(0)->±oo). (110) 
Proof. From the proof of Lemma A. 3 above, we know that 

" (n) W E f H ~t ') ' I™ e- ex x l+1 ^\x)Ax (111) 

for n=l,...,m. Using the assumptions (108) and the well-known relation between 
the asymptotic behavior of the integral modulus of continuity at zero and the asymp- 
totic growth of the Fourier transform at infinity yields (109) for n= 1, ...,m. To 
be more specific, we apply Kawata (1972), Theorem 2.7.4 to the functions f(x) = 
Q-^( 9 ) x x t + 1 v f t )(x)I( x> Qy Note that there is a misprint (not relevant here) in the ref- 
erence: f(t) should be f(t). 

The case n = follows immediately by using the estimate for k'(0) in n(6) = J Q k'(C) d£. 
Plugging these estimates into (87) and looking at the explicit formula for the Bell poly- 
nomials given in Comtet (1970), Theorem III. 3. A shows that 

PmiiP) = o(j2\ \ ai{a ~ 1)+ '" +am{a ~ m) ) =C(|6f ""'"). (112) 

□ 



Lemma A. 8. Suppose the assumptions for Lemma A. 7 hold and 

\(™)(0. t )= e K{e)t . (113) 

n dQ m dt n v ' 

Then 

\l™\8;t)=O(\0\ ( - rn+ ^ a - m ) (3(0)->±oo). (114) 
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Proof. First, we have 

X^(6;t)=K(e) n e K ^ t . (115) 

Using (89) and (109) with n = shows the claim for m = 0. Next, by differentiating (115) 
m > 1 times according to Faa di Bruno's formula, we get 

m 1 f 9\ 

^n l] (0;t) = X]E ( ) (n)M6) n ~ j *~ j e" m M6)<0) n e <(0t . (116) 
1=1 j=0 ^ ' 

Using (89), (109) with n = and (110), we obtain (114). □ 
A. 2. Convolutions and Laplace transforms 

In this subsection, we provide further auxiliary results for the proof of Theorem 3. For 
notational convenience, let us define V(x) = U + (x) and V n (x) = V* n (x). First it is shown 
in Lemma A. 9 that the convolution powers V* n (x) exist. We then show that the V* n (x) 
are n-times diffcrcntiable and we provide an integral representation in Lemma A. 10. 

Lemma A. 9. Let 

/>oo 

V(x)= u(y)dy. (117) 

J X 

Then 

V n (x)=V* n (x), x>0, (118) 
is well defined for n > 1 and we have the Laplace transforms 

' K((f\ n 

c- e ^„(x)dx-(-ir-yL, 5R(0)>O. (119) 
o & 

Proof. Let r > be arbitrary. A standard argument using the Fubini-Tonelli theorem 
shows that V(x) = e~ rx V(x) is intcgrablc. Thus, the convolution powers V* n (x) exist 
for almost all x > and are intcgrable on M>o. As we have 

e r*v* 2 (x)= [ V(y)V(x~y)dy, (120) 
Jo 

this shows that V(y)V(x — y) is intcgrablc on (0, x) and V* 2 (x) exists for almost all x > 0. 
Repeating the argument shows that the higher convolution powers exist for almost all 
x>0. □ 

In Barndorff-Nielsen and Hubalek (2006), a second proof is given, using the theory of 
convolutions of functions of the class Jo (from Doetsch (1950)) and it is shown that the 
convolution powers V n (x) = V* n (x) actually exist for all x > 0. 
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Next, we address the differentiability of V n (x) = V* n (x). For technical reasons, we 
work with the derivatives of x m V n (x) instead, where m is sufficiently large. 

Lemma A. 10. Suppose n € N and assumption A m {c) holds for some m > n + 2 and 
c>0. Then V n (x) is n-times differ entiable and 

(_l\m+n rc+ico / ..(a\n\ ( m ) 

V(*)( a .) = ijy / e e *d9, x>0, k = Q,l,...,n. (121) 

n y ' 2mx m J c _ ioo \e n - k ) y ' 

Proof. The result follows basically from Erdelyi et al. (1954), IV.4.1 (13), page 130, 
which implies, that 

e-^e)(x)d^=(^g) (m) . (122) 

For a rigorous proof of (121), we have to show that (i) V n (x) is actually n-times diffcr- 
cntiable, (ii) that the intcgrability conditions that allow the application of the quoted 
rule for the Laplace transform are satisfied and (iii) that the use of the Laplace inversion 
formula to deduce (121) from (122) is valid. This can be done in an elementary way by 
induction, with some careful and tedious bookkeeping, incorporating the estimates from 
Appendix A.l. For details, we refer the reader to Barndorff-Nielsen and Hubalck (2006), 
Lemmata 23 and 24. □ 



The result is less obvious than it at first seems. First, the common belief, that convo- 
lution increases smoothness is, in general, not true, as the shocking counterexamples in 
Uludag (1998) demonstrate. Second, the integrability assumptions to apply the standard 
theorems on the derivatives of convolutions (such as in Doctsch (1950), 1.2.14.5) are typ- 
ically not satisfied in our setting for V{x), as can be immediately seen from the positive 
stable example. Third, the application of Erdelyi et al. (1954), IV.4.1 (13), page 130 
cannot validly be decomposed into an application of Erdelyi et al. (1954), IV.4.1 (8), 
page 129 followed by an application of Erdelyi et al. (1954), IV.4.1 (6), page 129. 
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